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Measurement of interaction between antiprotons 


One of the primary goals of nuclear physics is to understand the force between nucleons, 
which is a necessary step for understanding the structure of nuclei and how nuclei inter¬ 
act with each other. Rutherford discovered the atomic nucleus in 1911, and the large body of 
knowledge about the nuclear force since acquired was derived from studies made on nucleons 
or nuclei. Although antinuclei up to antihelium-4 have been discovered ® and their masses 
measured, we have no direct knowledge of the nuclear force between antinucleons. Here, 
we study antiproton pair correlations among data taken by the STAR experiment ^ at the 
Relativistic Heavy Ion Collider (RHIC) ^ and show that the force between two antiprotons 
is attractive. In addition, we report two key parameters that characterize the correspond¬ 
ing strong interaction: namely, the scattering length (/o) and effective range (do). As direct 
information on the interaction between two antiprotons, one of the simplest systems of antin¬ 
ucleons, our result provides a fundamental ingredient for understanding the structure of 
more complex antinuclei and their properties. 

Although the theory of Quantum Chromodynamies (QCD) provides us with an understand¬ 
ing of the foundation of the nuclear force, this binding interaction in nuclei operates at low energy, 
where the force is strong and difficult to calculate directly from the theory. For that reason, a 
common parameterization of the effective interaction between nucleons is based on experimental 
measurements, and the corresponding parameterization for antinucleons remains undetermined. 
The important parameters in such a description of the interaction are the scattering length (/o), 
which is related to elastic cross sections, and the effective range of the interaction (do), which is 
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determined to be a few femtometers (the typieal nuclear scale). The existence and production rates 
of antinuclei offer indirect information about interactions between antinucleons, and also have rele¬ 
vance to the unexplained baryon asymmetry in the universe®. Antinuclei produced to date include 
antiprotons, antideuterons, antitritons, antihelium-3, and the recently discovered antihypertriton 
and antihelium-4 (see ref® and references therein). The interaction between two antinucleons 
is the basic interaction that binds the antinucleons into antinuclei, and this has not been directly 
measured previously. Of equal importance, one aspect of the current measurement is a test of 
matter-antimatter symmetry, more formally known as CPT — a fundamental symmetry of physi¬ 
cal laws under the simultaneous transformations of charge conjugation (C), parity transformation 
(P) and time reversal (T). While various prior CPT tests ® have been many orders of magnitude 
more precise than what is reported here, there is value in independently verifying each distinct 
prediction of CPT symmetry®. 

Ultra-relativistic nuclear collisions produce an energy density similar to that of the Universe 
microseconds after the Big Bang, and the high energy density creates a favourable environment for 
antimatter production. The abundantly produced antiprotons provide the opportunity to measure, 
for the first time, the parameters /o and do of the strong nuclear force between antinucleons rather 
than nucleons. 

The technique used to probe the antiproton-antiproton interaction involves momentum cor¬ 
relations, and it resembles the space-time correlation technique used in Hanbury-Brown and Twiss 
(HBT) intensity interferometry. Since its invention for use in astronomy by Robert Hanbury- 
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Brown and Richard Q. Twiss in the 1950’sl^, the HBT technique has been adopted in many areas 
of physios, inoluding the study of the quantum state of Bose-Einstein oondensates the oorrela- 
tion among eleotrons ^ and among atoms in oold Fermi gases ® A Bose-Einstein enhanoement in 
partiole physios was first observed in the late 1950’s as an enhaneed number of pairs of identioal 
pions produced with small opening angles (GGEP effeot) Eater on, Kopylov and Podgoretsky 
noted the common quantum statistics origin of HBT and GGEP effects, and, through a series of 
papers they devised the basios of the momentum oorrelation interferometry teohnique. In this 
teohnique, they introdueed the oorrelation funetions (CF’s) as ratios of the momentum distributions 
of eorrelated and unoorrelated partieles, P 2 ) = G = 1 for no oorrelations, 

suggested the so-oalled mixing teohnique to eonstruot the unoorrelated distribution by using parti¬ 
eles from different oollisions (events), and formulated a simple relation of the CF’s with the spaee- 
time struoture of the partiole emission region. As a result, the momentum oorrelation teohnique 
has been widely embraoed by the nuolear physios oommunity 
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Figure 1: A schematic of the two-particle correlation process in a heavy-ion collision overlaid on 
an event display from the Time Projection Chamber in the STAR detector. The curves show par¬ 
ticle trajectories, from which the track momenta are determined. These trajectories are measured 
in three dimensions, but are projected onto a single plane in this illustration. The STAR detector 
measures 3-vector momenta over a wide range beginning at about 0.1 GeV/c. Two particles emit¬ 
ted from two separated points, with four-coordinates Xa and X^, are detected with four-momenta 
Pi and p 2 . For the pair of indistinguishable particles with even/odd total spin, the two quantum me¬ 
chanical amplitudes (representing, for non-interacting particles, products of plane waves {pi\Xa) 
(p 2 \Xb) and (p 2 |Xa) (pi|Xb), where {p\X) = exp(—ipX)) must be added/subtracted to yield the 
amplitude which is symmetric/antisymmetric with respect to the interchange of particle momenta. 
This results in an enhancement/suppression in the joint detection probability at zero momentum 
separation with the width inversely proportional to the space-time separation of particle emission 
points. 
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Figure illustrates the process of constructing two-particle correlations in heavy-ion colli¬ 
sions. In addition to quantum statistics effects, final state interactions (FSI) play an important role 
in the formation of correlations between particles. FSI include, but is not limited to, the forma¬ 
tion of resonances, the Coulomb repulsion effect, and the nuclear interactions between two parti- 
ctes [Unmilll] fact, FSI effects provide valuable additional information. They allows for (see 
ref |i^|2Ql and references therein) coalescence femtoscopy, correlation femtoscopy with non-identical 
particles, including access to the relative space-time production asymmetries, and a measurement 
of the strong interaction between specific particles. The latter measurement is often difficult to 
access by other means and is the focus of this paper (for recent studies see refl^IEll). 

The data used here consists of Au -i- Au collisions at a centre-of-mass energy of 200 GeV 
per nucleon pair, taken during RHIC’s operation in the year 2011. In total, 500 million events 
were taken by the minimum-bias trigger at STAR. The minimum-bias trigger selects all particle- 
producing collisions, regardless of the extent of overlap of the incident nuclei, but with a require¬ 
ment that collisions must have occurred along the trajectory of the colliding Au ion and within 
±30 cm of the center of STAR’S Time Projection Chamber (TPC) 1^. Events are categorized by 
their centrality, based on the observed number of tracks emitted from each collision. Zero percent 
centrality corresponds to exactly head-on collisions which produce the most tracks, while 100% 
centrality corresponds to barely glancing collisions, which produce the fewest tracks. Events used 
in this analysis correspond to the 30%-80% centrality class, for which the signal due to two-particle 
interaction is the clearest. 
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The two main detectors used in the measurement are the STAR TPC and the Time of Flight 
Barrel (TOF) The TPC is situated in a solenoidal magnetic field (0.5 T), and it provides a 
three-dimensional image of the ionization trails left along the path of charged particles. The TOF 
encloses the curved surface of the cylindrical TPC. In conjunction with the momentum measured 
via the track curvature in TPC, particle identification (PID) is achieved by two key measurements: 
the mean energy loss per unit track length, (dE/dx), which can be used to distinguish particles 
with different masses or charges, and the time of flight of particles reaching the TOF detector. 
Figurej^shows a typical calculated mass-squared (m^) distribution, derived from the time-of-flight 
and tracking information, versus (see caption) for antiprotons. 
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Figure 2: m? versus for negatively charged particles. Here m? = (p^/c^)(f^c^/L^ — 1), 
where t and L are the time of flight and path length, respectively, c is the light velocity. 2 ; = 
\n{{dE/dx)/{dE/dx) e) and {dE/dx)E is the expected value of (dE/dx) for (anti)protons. az 
is the r.m.s width of the 2 ; distribution, and is the number of standard deviations from zero, 
the expected value of 2 : for (anti)protons. The antiprotons, centered at = 0.88 (GeV/c^)^ 
and = 0, are well separated from other particle species for the momentum range specified. 
(Anti)protons satisfying 0.8 (GeV/c^)^ < mass^ < 1 (GeV/c^)^ and < 1.5 are selected for 
making pairs. With this selection, the purity is > 99% for (anti)protons with transverse momentum 
less than 2 GeV/c. 

The population distribution of (anti)proton pairs as a function of (anti)proton momentum 
(k*) in the pair rest frame is measured for the correlated pairs from within the same event, A{k*), 
and, separately, for the non-correlated pairs from two different (mixed) events, B{k*). The ratio of 
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the two, gives the measured CF (see Methods). The observed (anti)protons ean eome from 
weak deeays of already eorrelated primary partieles, henee introdueing residual eorrelations whieh 
eontaminate the CF. The dominant eontaminations to the CF eome from the p-A (p-A) and A-A 
(A-A) eorrelations, and are taken into aeeount by fitting the CF with eorresponding eontributions. 
Taking the two-proton eorrelation measurement as an example^, 

CindnsWeik*) = 1 + Xpp[Cpp{k* ] Rpp) - 1] + XpA [CpA ^ “ 1] + ^AA[CAA(fc*) “ 1], (1) 

where Cindusiveik*) is the inelusive CF, and Cpp{k*; Rpp) is the true proton-proton CF, whieh ean 
be deseribed by the Lednieky and Lyuboshitz analytieal model In this model, for given s- 
wave seattering parameters, the eorrelation funetion with FSI is ealeulated as the square of the 
properly symmetrized wave funetion averaged over the total pair spin and the distribution of rel¬ 
ative distanees of partiele emission points in the pair rest frame (see Methods). Cp\{k*] Rp^) is 
the proton-A CF from a theoretieal ealeulation whieh ineludes all final-state interaetions and 
explains experimental data well^. Cw^k*) is from an experimental measurement eorreeted for 
misidentified A’s^^. Rpp and Rp\, assumed numerieally to be the same, are the invariant Gaussian 
radii from the proton-proton eorrelation and the proton-A eorrelation, respeetively. Xpp, Xp\ and 
xaa, taken from the THERMINATOR2 model are the relative eontributions from pairs with 
both daughters from the primary eollision, pairs with one daughter from the primary eollision and 
the other one from a A deeay, and pairs with both daughters from a A deeay, respeetively. 
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Figure 3: Correlation function for proton-proton pairs (top), antiproton-antiproton pairs (middle), 
and the ratio of the former to the latter (bottom). Errors are statistical only. The fits to the data with 
equation!^ C-mciusiveik*), are plotted as solid lines, and the term l + Xpp[Cpp{k*; Rpp) — 1] is shown 
as dashed lines. The x^/ndf of the fit is 1.66 (1.61) for the top (middle) panel. To take advantage 
of the existing knowledge on the proton-proton interaction, which is relatively well understood, 
when fitting the proton-proton correlation, /o and do for protons are fixed at values measured from 
proton-proton elastic-scattering experiments, which are 7.82 fm and 2.78 fm, respectively^. When 
fitting the antiproton-antiproton correlation, /o and do are treated as free parameters. 


Figure shows the CF for proton-proton pairs (top panel) and antiproton-antiproton pairs 
(middle panel), for the 30 — 80% centrality class of Au -i- Au collisions at a centre-of-mass energy 
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of 200 GeV per nucleon pair. The proton-proton CF exhibits a maximum ai k* ~ 0.02 GeV/c 
due to the attractive singlet s-wave interaction between the two detected protons and is consistent 
with previous measurements 1^. The antiproton-antiproton CF shows a similar structure with the 
maximum appearing at the same k* value. In the bottom panel, the ratio of the inclusive CF for 
proton-proton pairs to that of antiproton-antiproton pairs is presented. It is well centered at unity 
for almost all the k* range, except for the region k* < 0.02 GeV/c, where the error becomes 
large. This indicates that the strong interaction is indistinguishable within errors between proton- 
proton pairs and antiproton-antiproton pairs. By fitting the CF with equation we determine the 
singlet s-wave scattering length and effective range for the antiproton-antiproton interaction to be 
/o = 7.41±0.19(stat)±0.36(sys) fmanddo = 2.14±0.27(stat)±1.34(sys) fm, respectively. The 
extracted radii for protons {Rpp) and that for antiprotons {Rpp), are 2.75 ± O.Ol(stat) ± 0.04(sys) 
fm and 2.80 ± 0.02(stat) ± 0.03(sys) fm, respectively. 
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Figure 4: The singlet s-wave scattering length (/q) and the effective range (do) for the antiproton- 
antiproton interaction plotted together with the s-wave scattering parameters for other nucleon- 
nucleon interactions. Here, statistical errors are represented by error bars, while the horizontal 
uncertainty for /q is smaller than the symbol size, and systematic errors are represented by the 
box. Errors on other measurements I2ZI22! are on the order of a few percent, smaller than the symbol 
size. 


Figure presents the first measurement of the antiproton-antiproton interaction, together 
with prior measurements for nucleon-nucleon interactions. Within errors, the /o and do for the 
antiproton-antiproton interaction are consistent with their antiparticle counterparts - the ones for 
the proton-proton interaction. Our measurements provide parameterization input for describing the 
interaction among cold-trapped gases of antimatter ions, as in an ultra-cold environment, where s- 
wave scattering dominates and effective-range theory shows that the scattering length and effective 
range are parameters that suffice to describe elastic collisions®. The result provides a quantitative 
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verification of matter-antimatter symmetry in the important and ubiquitous context of the forces 
responsible for the binding of (anti)nuclei. 

Methods 

Event mixing for non-correlated pairs and the correction for purity. Non-correlated pairs each 
consist of two daughters particles. These daughters belong to two events which are carefully 
chosen so that they have similar event multiplicity and topology. The ratio of the (see above), 
after being normalized at a large k* (at least 0.25 GeV/c), gives the measured CF, G(fc*)meas- 
Because in practice one cannot select 100% pure (anti)protons, a correction to pairs is applied to 
obtain the PID-purity-corrected CF: GpurityCorrected(fc*) = simplicity, in Eq.[^ 

the subscript “meas” is dropped, and elsewhere in this paper, the subscript “PurityCorrected” is 
dropped. 

The transformation from k*j^ to k'^. Cpx{k*) in equationis naturally expressed as a function 
of k*j^. Thus, to use it in the above equation, one needs to transform it to a function of k*p. Here 
k*j^ (and k^) is the magnitude of the three-momentum of either particle in the pair rest frame, 
while in this case for k*p, one of the protons is the decay daughter of A. This transformation is 
done by CpA{k*pp) = jCpA{k*j^)T{k*j^, k*p)dk*^, where T{k*^, k*p) is a matrix that transforms 
k*A to The transformation matrix is generated with the THERMINATOR2 model ^ which 

is a Monte Carlo event generator dedicated to studies of the statistical production of particles in 
relativistic heavy-ion collisions. 
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The calculation of the FSI contribution to the correlation function. The femtoscopic corre¬ 
lations due to the Coulomb FSI between the emitted electron and the residual nucleus in beta 
decay have been well known for more than 80 years; they reveal themselves in a sensitivity of the 
Fermi function (an analogue of the CF^^ to the nuclear radius. Compared with non-interacting 
particles, the FSI effect in a two-particle system with total spin S manifests itself in the sub¬ 
stitution of the product of plane waves, exp(—— ip 2 Xh), by the non-symmetrized Bethe- 
Salpeter amplitudes ^b) = X^) | i^|i9|32 | 33 | p^j. identical particles, the sym- 

metrization requirement in the representation of total pair spin S takes on the same form for both 
bosons and fermions: the non-symmetrized amplitude should be substituted by X;,) -f 

{-lY^^^-}{Xa,Xi,)]/y/2. In the pair rest frame, Xa-Xb = {C, r*} andpi-p 2 = -Wa) 2k*} 
where oj* = {mf + is the energy of a particle of mass mj, and t* and r* are the rela¬ 

tive emission time and relative separation in the pair rest frame, respectively. In this frame, the 
non-symmetrized Bethe-Salpeter amplitude at equal emission times (t* = 0) reduces, up to an 
inessential phase factor, to a stationary solution of the scattering problem, At small 

relative momenta, k* < ~l/r*, this solution can be used in practical calculations with the con¬ 
dition |f*| ^ The equal-time approximation is almost exact in beta decay, and it is 

usually quite accurate for particles produced in high-energy collisions (to a few percent in the FSI 
contribution to CF’s of particles even as light as pions^^. In collisions involving heavy nuclei, the 
characteristic separation of the emission points, r*, can be considered substantially larger than the 
range of the strong-interaction potential. The FSI contribution is then independent of the actual 
potential form and can be calculated analytically with the help of corresponding scattering ampli- 
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tudes only® At small k*, it is basically determined by the s-wave scattering amplitudes f^{k*) 
sealed by the separation r* 

The analytical calculation of the (anti)proton-(anti)proton correlation function. The (anti)proton- 
(anti)proton eorrelation function, Cpp{k*; Rpp) in equation[^ can be deseribed by the Lednicky and 
Lyuboshitz analytieal model®. In this model, the correlation funetion is calculated as the square of 
the properly symmetrized wave funetion averaged over the total pair spin S and the distribution of 
relative distances (r*) of particle emission points in the pair rest frame, assuming 1/4 of the singlet 
and 3/4 of triplet states and a simple Gaussian distribution dN/d?Y* ~ exp(— r*^/4i?pp). Start¬ 
ing with the FSI weight of nueleons emitted with the separation r* and deteeted with the relative 
momentum k*, 

m(k*,r-) = |^!it'((r-) + (-l)*V-/!+>(r*)|V2, 

where (r*) is the equal-time {t* = 0) redueed Bethe-Salpeter amplitude whieh ean be ap¬ 
proximated by the outer solution of the seattering problem®. This is 

where rj = {k*ac)~^, ac = (57.5 fm) is the Bohr radius for two protons, p = k*r*, ^ = k*r* -f p, 
Ac{r]) is the Coulomb penetration faetor given by Ac{r]) = 2Tir][exTp(2Tir]) — 1 ]“^, F is the eonfluent 
hypergeometrie funetion, G{p, rj) = ^JaJ^[Gq{p, rj) + iFo^p, rj)] is a eombination of the regular 
(Fo) and singular (Go) s-wave Coulomb funetions, 

fcik*) = [^ + Tydok*"^ - —h{T]) - ik*A^{T])]-^ 

Jo ^ 
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is the s-wave scattering amplitude renormalized by the Coulomb interaction, and h{r]) = rf Y, 

n=l 

— C — In \r]\ (here C = 0.5772 is the Euler constant). The dependence of the scattering pa¬ 
rameters on the total pair spin S is omitted since only the singlet (S = 0) s-wave ESI contributes 
in the case of identical nucleons. The theoretical CE at a given k* can be calculated as the average 
ESI weight (r(;(k*, r*)) obtained from the separation r*, simulated according to the Gaussian law, 
and the angle between the vectors k* and r*, simulated according to a uniform cosine distribution. 
This CE is subject to the integral corrections^ —y4c(?7) |/c(fc*) pdo/(8-\/7r-Rpp) due to the deviation 
of the outer solution from the true wave function in the inner potential region. In addition, in 
Au-i-Au collisions the emitting source has a net positive charge, and it influences the CE differently 
for proton and antiproton pairs. This effect is included in the consideration according to ref. 

Systematic uncertainties. The systematic uncertainties include variations due to track-wise and 
pair-wise cuts, the uncertainty in describing the correlation function SH, and the uncertainty 
from the Caa measurement. The latter dominates the systematic error of do and /o, and it affects 
do more than it does /o because the shape of the CE is sensitive to do, in particular at low k*. As 
a consistency check, when fitting the proton-proton CE, both /o and do are also allowed to vary 
freely, and the fitted /o and do agree with the results from fitting the antiproton-antiproton CE. 
Assuming the measurements from different systematic checks follow a uniform distribution, the 
final systematic error is given by (maximum - minimum)/A/l2. In our calculations, we consider 
the two-proton wave function, taking into account the Coulomb interaction between point-like 
protons in all orbital angular momentum waves and the strong interaction in the s-wave only. 
We neglect the small non-Coulomb electromagnetic contributions due to magnetic interactions. 
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vacuum polarization, and the finite proton size HZOZlIlIl. This approximation ehanges the seattering 


parameters at the level of a few pereentl^ZHUlll. 
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